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ABSTRACT 

Fibre modal noise occurs in high spectral resolution, high signal-to-noise applications. 
It imposes fundamental limits on the photometric accuracy of state-of-the-art fibre- 
spectrograph systems. In order to maximize the performance of current and future 
instruments it is therefore essential to predict fibre modal noise. To attain a predictive 
model we are using a dual approach, bringing theoretical assumptions in line with the 
experimental data obtained using a test-bench spectrograph. We show that the task 
of noise prediction can be reduced to determining the visibility of the modal pattern 
which can be measured at the detector plane. Subsequently, the visibility-dependence 
of essential parameters is presented. This work will soon provide a basis for prediction 
of modal noise limitations in fibre-coupled spectrograph designs. 

Key words: instrumentation: spectrographs, techniques: spectroscopic, methods: 
laboratory, methods: analytical, line: profiles, techniques: radial velocities 



1 INTRODUCTION 
1.1 Modal noise 

Modal noise is a measurement uncertainty caused by in- 
herent statistical properties of an optical fibre. It can be 
ascribed to an interference effect at the fibre output, oc- 
curring when light of limited bandwidth is detected. Hence, 
modal noise is particularly important for high-resolution ap- 
plications where central wavelength determination, spectral 
line ratios and shapes play a crucial role. This paper concen- 
trates on the effect of modal noise on the photometric ac- 
curacy in a measurement. Generally speaking, any mode fil- 
tering process affects the throughput of the combined fibre- 
spectrograph system. In particular, beam truncation (com- 
pare fig. [TJ will limit the achievable signal-to-noise of the 
instrument. 

As telescope apertures increase in size and ever higher spec- 
tral resolving powers become available, higher signal-to- 
noise ratios (SNR) can be achieved and the light at the de- 
tector becomes more coherent. The employment of optical 
fibres plays an important role here in order to achieve the 
required instrument stability and spectral resolving power 
(see |Heacox| ( |l988| ) and |Allington-Smith| pOOo) ). Unlike in 
conventional slit spectroscopy, fibres only partially retain 
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spatial and angular information. For differing seeing condi- 
tions the light that is coupled into the fibre largely causes the 
same spectrograph response and thus stabilizes the signal; 
A property that is referred to as scrambling in the literature 
(see Avila & Singh (2008J). With these advantages and the 
increasing requirements on instrument precision, the modal 
noise issue will continue to gain in prominence. It is there- 
fore essential to predict modal noise in order to determine 
the performance of the complete instrument, allowing for 
optimization and precautionary measures that help reduc- 
ing modal noise before it becomes a fundamental constraint 
to the instrument. 



Epworth ( 1979 1 first noted that the modal properties of 



fibres will limited the accuracy of intensity measurement. 
Inherent fibre properties cause a speckle pattern at the fi- 
bre exit which is subject to changes in the fibre condi- 
tion, e.g. stress or temperature shifts. In a subsequent paper 



by Rawson et al. ( 1980p the bandwidth-dependence of the 
speckle contrast was studied: A lower coherence (e.g. in low- 
resolution spectrographs) results in diminished speckle con- 
trast (low visibility). This theory predicts a strong length- 
dependence which was disproven for fibre lengths < 10 2 m 
by |Baudrand fc Walker| ( p00l| ; |Corbett| pOOo] ). Although 
coherence effects are a key parameter, attention has to be 
paid to the statistical properties of the speckles and their ef- 
fect on modal noise. As was recently pointed out by |Grup"p| 
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(20031 spatially filtering apertures are often deployed, e.g. 



fibre-slit at spectrograph entrance or overfilled grating, ini- 
tially aimed at high spectral resolution at low grating cost - 
but with disregard to limiting modal noise. 
The first investigation to address aperture clipping was con- 
ducted by |Hill et al. (|1980[), and soon followed by Good- 



man & Rawson ( 1981 1. This work showed that modal noise 



is caused by mode filtering processes but is restricted to 
the extreme case of highly coherent light. The first work to 
combine coherence effects and speckle statistics into a con- 
cise model is reported by Kanada & Aoyama (119831). Their 



experimental results show good agreement with theoretical 
predictions, but the investigation concentrates on high con- 
nection losses, i.e. strong beam vignetting. Moreover, the 
considered visibility regime of > 10% is significantly higher 
than relevant in astronomical spectroscopy (5% and less in 
this work) . The first prediction of the wavelength-dependent 
spectrograph response was achieved by |Chen et al.| ( 2006 ) 
for few-mode fibres. Their mode-propagation simulation is 
in excellent agreement with experimental findings, allowing 
for precise modal noise prediction. But the multi-mode case 
is more complicated and already the absence of strong fibre- 
length dependence suggests that our theoretical picture of 
light propagation in optical fibres is not yet sufficiently de- 
veloped. 

Today's instruments are already meeting challenging re- 
quirements, among these the most outstanding are instru- 
ments dedicated for exoplanet research using the radial ve- 
locity method. High-resolution systems like HARPS and 
ESPRESSO are now striving to detect Earth-like planets, 
but the unprecedented precision required here is found to 
be restricted (amongst other effects) by modal noise (see 
|Wilken et ah] p010| ; poisse et al.| poTo] )). 
Stellar observations that measure metallicities or make use 
of analytical tools like the Stark effect for temperature de- 
termination will be affected. Also studies on interstellar ma- 
terial measuring ionization levels will ultimately be limited 
by modal noise. Furthermore, for solar observations working 
at high resolution, high signal-to-noise levels are affected by 
modal noise (e.g. high-precision-polarimetry) . 



1.2 Scope of this investigation 

This paper addresses modal noise in the parameter range 
of spectrographs dedicated to astronomical applications. We 
concentrate on intermediate to low aperture restrictions and 
visibilities < 10%. In contrast to the literature discussed 
above, we consider the effect of modal noise for the low co- 
herence case (in contrast to the investigations at high coher- 
ence found in literature) and for the fibre spectrum rather 
than just the fibre image. 

A full theoretical picture of the effects involved, requires 
understanding phenomena like intermodal coupling, modal 
filtering and spatial mode distribution. Therefore, we feel it 
necessary that models with the objective of modal noise pre- 
diction need to be grounded on experimental results first. 
We start out in section [2] with the existing theoretical back- 
ground of modal noise. A first simplification is possible with 
the assertion (which will be experimentally verified), that 
modal noise can be described as the product of two indepen- 




P 2 = A d /A f 



Figure 1. Modal noise originates from truncating the beam of 
an inhomogeneous, timely varying intensity distribution. The es- 
sential geometrical variable p 2 is defined by two parameters; Af 
is the area of the fibre image, whereas is the part of the illu- 
minated area that is detected. 



dent parameters; a visibility and a coherent term. Hence, 
it can be separated into two sub-problems: The coherent 
speckle term is well understood from fibre and restricting 
aperture geometry, whereas visibility is mainly dependent on 
resolving power but also on the slit width. A simple model 
for visibility is derived from simplified ray propagation as- 
sumptions, outlining basic dependencies on spectrograph pa- 
rameters. 

For the experimental investigations a medium to high reso- 
lution (R = A/ A A = 10 000 - 200 000) Echelle spectrograph 
test-bench was built. Its design and the data-analysis are 
presented in section [3] The experimental findings (section 
4.1 1 support the above hypothesis about the modal noise 
factorization, the coherent part is proven to follow the theo- 



retical predictions of Goodman & Rawson (1981). The prob- 



lem of modal noise prediction thus reduces to the problem of 
determining the visibility (section |4.2[ | . This fact can already 
be used to directly extract modal noise information from 
raw data of calibration sources. To enable direct prediction 
of visibility and modal noise, we furthermore investigate the 
dependence of visibility on spectrograph parameters such as 
resolving power and input focal ratio. Although some agree- 
ment between theoretical predictions and experimental data 
was shown, additional investigation needs to follow to refine 
the visibility-model and thus enable modal noise prediction. 
This work will be subject to a later paper (Lemke et al., in 
preparation). 



2 THEORETICAL PREDICTIONS 

In confined waveguides as represented here by the standard 
optical fibre, the electromagnetic field has to fulfill certain 
boundary conditions. As a result, only the fraction of 
the light that matches the potential walls can propagate 
through the fibre and is quantized into fibre modes. The 
typical fibre differs from the ideal circular symmetric 
waveguide - along the length of the fibre a variation of 
parameters like core-diameter and refractive index occurs. 
These statistical properties cause a change in symmetry 
and phase relation of the supported modes, a phenomenon 
that manifests in a speckle pattern as in fig. [T] visible in 
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should serve us as a preliminary theoretical framework. For 



coherent light. Changing stress in the fibre (as for instance 
due to telescope tracking) causes the speckle pattern to 
change and in connection with mode filtering processes 
such as beam truncation causes modal noise. 
The SNR derives from theoretical considerations as follows: 
According to |Goodman fc Rawson| ( |1981[ ) the SNR can be 
calculated from speckle statistics and geometrical consider- 
ations, given the number of speckles (which is equal to the 
number of excited modes M, see |Daino et al.| (1980)) and 
the area fraction p 2 of the illuminated fibre end-face that is 
transmitted to the detector surface (see fig. [I]): 



SNRoah 



M + l 



(1) 



Note that the truncation described by p 2 is not necessarily 
generated by a slit at e.g. the spectrograph entrance, but 
also by overfilling the grating or any other spatially filtering 
process (e.g. inhomogeneous sensitivity of the CCD). The 
subscript coh indicates that we consider the coherent case. 
Note that the question of how many modes are excited is 
not trivial. Although the mode volume V for a given input 



illumination can be directly calculated (see e.g. Hill et al. 
( [1980] )) from the fibre core radius a, wave number k = 2tt/X 
and angle 6i n of the light cone feeding into the fibre 

V = i (ak6 m ) 2 , (2) 

this is only an upper limit. M very much depends on the 
fibre-launch geometry (axis-offset, angle), as was asserted 
by |Corbett eFaL] ( [2007] ) . 

However, these considerations only hold true when consid- 
ering coherent light. The speckle contrast is reduced when 
observing light of broad bandwidth and studying a large set 
of fibre modes, as is the case for fast coupled multimode 
fibres. Then the fibre is populated by a large number of 
modes, resulting in a broad range of propagation constants. 



This is taken into account by Kanada & Aoyama ( 1983 1, di- 



viding the pure speckle statistics expression by the visibility 
v, yielding 

SNR=-SNR coh . (3) 
v 

Here v can be theoretically derived by taking into ac- 
count the finite bandwidth of the optical source (/ c ) and the 
phase delay between modes (~/ s ): 

-1/2 



t'th 



= [1 + 2(f s /f c f] ■ 



(4) 



An approximation (see Appendix A) helps understand- 
ing the structure of the equation and its dependence on 
spectrograph parameters, leading to the following result: 



Vth = 



1 + 2 



L6 2 
2n c RX c 



-1/2 



(5) 



for fibre length L, numerical aperture 9, core refractive-index 
n c , spectrograph resolving power R and central wavelength 
Ac 

We make no claim here that these simple assumptions lead 
to a model fully consistent with the experiment, but they 



instance it was shown by Baudrand & Walker (20011; Cor 



bett (20061, that the SNR has no significant dependence 



on fibre length for L 100 m. This is possibly because the 
fibre does not terminate with a perfect quartz-air intersec- 



tion. This so called end-effect (Poppett & Allington-Smith 



( |2010[ )) may lead to a stronger mode mixing and thus less co- 
herence between the modes, which would explain the dimin- 
ished visibility. In the following L will therefore be treated 
as a free parameter. Thus, v remains a function of 8, n c , R, 
X, we will concentrate here on 9 and R only. In the limit of 
low visibilities it applies that (L8 2 ) / (2n c RX c ) 2> 1, so that 
eq. (JsJ can be approximated by 

vth oc 6r 2 s-\ ( 6 ) 

where R is assumed to be proportional to the slit width 
s. We will see that this approximation is fulfilled for the 
#-case and fulfilled at least by trend for s, whereas the 
full dependence of visibility on the spectrograph resolving 
power R will remain a matter for further investigations. 



3 EXPERIMENTAL SETUP 
3.1 Spectrograph setup 

The experimental apparatus consists of fibre input, 
agitator, slit unit and spectrograph. 



fibre 



Fibre input (fig. |2(a)[ |; The test fibre is illuminated 
using a power stabilized halogen lamp behind a diffuser 
plate and a 500 pm pinhole. The lens system demagnifies 
the spot to a image size of 80 to 90 pm at the fibre input 
plane. An iris is used to stop the beam down in order 
to control the input focal ratio F. Focal ratios F = 2.6, 
3.9, 6.5 and 13.2 were used, by default F = 3.9 if not 
stated otherwise. The high focal ratio of 3.9 is to mimic 
astronomical devices where the fibre is fed with a fast beam 



to minimize focal ratio degradation (FRD, Ramsey ( 1988 1 ) , 
an effect causing the output focal ratio F ou t to be generally 
smaller than the input focal ratio. Note that the focal ratio 
notation (F or f/#) is applied here (and widely in optics 
literature) whereas often the geometrically more descriptive 
6 = arctan(2F) is used in order to present the functional 
dependencies more clearly. 

It is important here that the fibre core is fully illuminated. 
Choosing a significantly smaller spot-size results in poor 
scrambling. Incomplete scrambling manifests itself as a 
global structure of the intensity pattern in the spectrum 



(see Avila & Singh ( 2008 1). The pattern is dependent on 
the position of the input spot, the variation exclusively 
occurs in cross-dispersion direction and thus differs from 
modal noise. This inhomogeneous intensity distribution 
is dependent of the position of the illuminating input 
spot. We chose a larger spot-size (ss 70 % of the fibre-size) 
which resulted in a more homogeneous intensity pattern. 
A microscope is used to view the lateral position and the 
focus position of the pinhole image on the fibre core. 
The beamsplitter is specified as "broadband non-polarizing" 
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Figure 2. (a) Fibre input: The focus position can be controlled with the aid of a microscope in order to obtain optimal fibre illumination; 
(b) Fibre agitation scheme: The fibre is fixed at two points with a large free length in between. An agitator is used to produce large 
deflections in this region. 




Figure 3. Schematic depiction of the Echelle spectrograph used for experimental investigations. The fibre far-field is imaged onto the 
slit plane which coincides with the spectrograph entrance. High spectral resolution is achieved by using the Echelle grating in Blaze 
configuration. Light is filtered through a narrow bandpass filter to avoid spectra overlapping at the CCD. (Slit and 4 mm lens are not to 
scale, grating angle exaggerated.) 




(a) (b) 

Figure 4. (a) Re-imaged at the CCD is the fibre far-field as it occurs at the slit plane at the spectrograph entrance. Shown are the far- 
fields (slit removed) at different strengths of fibre-agitation (from left to right): No agitation; agitation 1.5 Hz; agitation 80 Hz; additional 
manual agitation. Exposure time was 3 s; (b) The fibre (far field) image can be truncated by a slit only as in fig. ^ or by additional 
apertures that truncate the beam in cross-dispersion direction, e.g. when restricting the slit-aperture accordingly (as in section [4.1[ l or 
when overfilling the dispersive grating. 
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Table l.Slit widths, Resolving power R, and pixel per resolution lem 5 th h = 800 mm (/i = 300 mm for lower resolution), 

element p for the two cases of different focal length. the fibre end-face is re-imaged onto the Echelle grating. 

The Echelle grating is used in Blaze configuration, with 
the reflection angle ip = 18°. A second lens of focal length 
ji = 800 mm (ji = 400 mm for lower resolution) images 
the slit-restricted fibre far field onto the CCD. The spectral 
dispersion has been measured using a Neon lamp and the 
width of one spectral resolution element was determined 
using laser illumination. With these measurements the 
spectral resolving power can be calculated (see table [I]). A 
fold mirror is applied in order to fit the spectrograph on 
the optical bench. Mirrors and gratings are known to have 
polarizing characteristics and thus potentially attenuating 
a portion of the modes. Therefore the polarization was 
measured by introducing a polarizing sheet after the slit 
plane (see fig. [3]l and found to be lower than 4%. Together 
with the beam splitter, we concluded that the influence 
of polarizing effects resulting from mirror and grating 
are marginal and only lead to a small deviation from the 
non-polarized case (appendix B). Therefore, the SNR in 
the following sections are the directly measured values, no 
polarization correction has been applied. Note that this 
needs to be taken into consideration when predicting modal 
noise effects in highly polarization filtering applications, 
e.g. when describing modal noise in polarimeters. 
In order to capture the spectra, an Atik Instruments 
314L+ camera was used. A narrow band filter (HeNe, 
632.8 ± 3.0 nm) is placed directly in front of the CCD 
to make sure that spectra of adjacent orders would not 
overlap. In addition, the slit-unit and spectrograph are 
placed in a light-tight box to keep the background signal 
low and air currents to a minimum. 



by manufacturer's specifications, 
the fibre launch, it was shown 

non-polarizing for the wavelength range under investigation. 
Utilizing a polarizing sheet before the beamsplitter, a po- 
larization signal could not be detected beyond the detector 
noise, the beamsplitter was thus found to be non-polarizing 
(less than 10% polarization for unpolarized input light). 

Fibre agitation and slit unit (fig. 2(b) and fig.[3|; After 
the fibre input, the fibre is attached to a loudspeaker allow- 
ing the fibre to be agitated (for flat-field image acquisition). 
The agitation procedure can be understood as applying 
temporarily varying stress that is applied to the fibre. 
The refractive index in the fibre thus changes accordingly, 
causing the speckle pattern at the output to vary. For 
integration times S> the period of agitation the recorded 
image is a pattern that (ideally) does no longer contain 
any modal information. These are used for flat-fielding (see 
section 3.2 1. 

In order to obtain the most efficient agitation, 



the fibre is 

secured on two points where the second point of securing is 
chosen at a larger distance from the point of agitation (see 
fig. [2JI. This allows the fibre to oscillate around its point 
of rest at this end. As it is believed that a large deflection 
changes the fibre properties by a larger amount, the distance 
is adjusted to the frequency at which it produces standing 
waves. An agitation frequency of 80 Hz (fig. |4(a)"| gives the 
best agitation results (low SNR in the respective exposures), 
and the modal pattern achieves a low contrast. Operating 
the agitation device at higher frequencies does not visibly 
improve the contrast, nor could a significant change in 
SNR be observed. From laser exposures it appears that 
additional manual agitation would improve the result, but 
again no difference for the SWR-value could be measured. 
The fibre end-face is placed behind a lens with focal length 
/o = 4 mm, imaging the fibre far-field on the slit plane 
where the slit holder can carry slits of the widths 50, 100, 
200 or 500 um, realizing the spectral resolving powers as 
documented in Table[T] Note that placing a slit at that 
plane causes beam vignetting so that p 2 < 1, thus inducing 
modal noise (cf. eq. [I]). 

Echelle spectrograph (fig.j3|; The spectrograph features 
a 4f-configuration, popular in astronomical applications 
(image and pupil planes are situated in the focal planes of 
the imaging lenses), this has the advantage of homogeneous 
illumination of the grating, minimizing aberrations and 
allowing for stray light filtering. Using a lens with focal 



3.2 Data analysis 

1. Modal noise determination 

Each dataset consists of a set of exposures recording 
the fibre at rest (typically six images) containing the 
inherent modal noise signature of the bare fibre. These are 
supplemented by a set of exposures recording the fibre in 
the agitated state. Additionally, background images of the 
same exposure time are recorded. For short exposure times 
the SNR reaches only low values, while longer exposure 
times result in increasing SNR approaching the theoretical 
limiting S/VR-level. Fig. [5] depicts how the exposure time 
was chosen: For increasing exposure times the data was 
recorded and its SNR determined until a plateau value was 
reached. This data point is then used for further analysis. 
Choosing even longer exposure times diminishes the SNR 
which indicates the occurrence of an additional source of 
noise. Very long exposure times can thus be critical and 
lead to a diminished SNR that is overestimating the modal 
noise. In fact, we assume here that the limiting SNR- value 
is equal to modal noise, although other sources of noise can 
contribute significantly, especially in the high SVR-regime. 
Potential sources of noise are Poisson noise (eliminated for 
long exposures), variations in the light source intensity and 
colour as well as fluctuations in the speckle pattern during 
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Figure 5. Exposure time was increased and for each dataset 
the modal noise calculated in order to find the point where the 
limiting signal-to-noise is reached. In this example 6 s was chosen 
for the exposure time of the image stack, exposure times increase 
for higher signal-to-noise ratios. 



collection of the unagitated fibre spectra, which could be 
caused by air currents (minimized) or shift in temperature. 
The error in the SNR of the unagitated fibre images of the 
measurement in fig. [5] is estimated during the data analysis 
process and of the order of 30-40, too small to appear in 
fig. [5] In general, SNR measurements are conducted at least 
two or three times and their variation is taken as a measure 
for the uncertainty of the measurement and analysis process. 

The data analysis involves basic reduction steps: 

• Background subtraction 

• Spectra extraction 

• Flat-fielding 

• Calculate standard deviation 

The background was subtracted using the correspond- 
ing dark-images. For exposure times longer than 10 s the 
images exhibit hot pixels ^ 3 times the average background 
signal; typically around five in the entire spectrum, despite 
the camera being cooled to 0°C. These pixels are retrieved 
in the dark-images and removed in all exposures by replac- 
ing each with the mean- value of adjacent pixels. 
We extracted the spectra by integrating over a large 
enough area in the cross-dispersion direction (300 pixels) 
for focal ratios up to 200 pixels full- width-half- maximum; 
the influence of background compared to the bright fibre- 
spectra is negligible. This could be shown from modal 
noise calculations using different sizes of integration ranges. 
Also, the read-out noise is <4ADU (i.e. 1 photo electron) 
per pixel (manufacturer's specifications) and is therefore 
neglected. 

After spectra extraction we are left with a one-dimensional 
array containing intensity information over spectral posi- 
tion. In order to remove any instrumental signature (e.g. 
beating pattern due to the bandpass-filter) , the unagitated 
fibre exposures are flat-fielded against the exposures of the 
agitated fibre. The quality of the flat-fielding is monitored 
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Figure 6. Unagitated fibre images are flat-fielded yielding a 
modal pattern. For visibility determination, adjacent minima and 
maxima are selected (circles) and the intensity of five pixel added 
(cross-shaped area) to feed into eq. |7| . The linear greyscale is 
normalized to the mean intensity. 



by determining the modal noise for an independent set of 
images of the agitated fibre; This is compared with the 
actual modal noise as recorded by the unagitated fibre 
exposure. 

The flat-fielded spectra are binned according to the size of 
one spectral resolution element (see Table [T]). In the final 
step the standard deviation between the bins is calculated 
in order to obtain the signal to noise value. Note that in the 
case where the number of modes is small (as presented by 



Chen et al. (20061), the likelihood of measuring a certain 
intensity does not resemble the usual error function. In 
this case a Bayesian approach is necessary. However, as 
we are concerned here with the highly multi-mode case, 
the central limit theorem can be applied and the intensity 
values scatter in the familiar Gaussian distribution around 
the central value; i.e. the standard deviation is a valid 
approximation for the one sigma confidence level and error 
propagation law can be applied as usual. 
After calculating the modal noise for the unagitated fibre 
exposures, the same procedure was applied to the exposures 
of the agitated fibre to compare with the theoretical photon 
noise, which was directly derived from Poisson statistics. 
With a few exceptions, all agitated fibre modal-noise values 
match with the theoretical photon noise, while always 
achieving significantly higher STVR-values than for the 
unagitated fibre. The remaining S/VR-difference to the 
Poisson value is probably due to additional sources of noise 
(as discussed above). 

2. Visibility determination 

The image in fig. [5] has been created by flat-fielding the 
central core area of the fibre-spectrum before applying the 
spectral extraction to the data. The pattern clearly shows 
a modal structure, the origin for modal noise. The visibility 
in each of these speckle patterns is derived by calculating 
the contrast between neighboring minima and maxima at 
six different positions of the spectrum: 



+ In 



(7) 
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Figure 7. Signal to noise and visibilitity vs slit sizes s. 
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Figure 8. Signal to noise and visibility vs input focal ratio < 



From these six visibility values the mean value is calculated, 
their scattering is used to estimate the uncertainty. 



4 EXPERIMENTAL RESULTS 

Fig.[7]and fig.[8]show slit size s and input focal ratio fl,„ vs. 
signal-to-noise SNR and visibility v, respectively. These val- 
ues are directly derived as described in section |3~2| For this 
particular setup, the 52VR-values are of the order of 10 3 and 
lower, and the visibility-values almost without exception do 
not exceed 5%, typical values are found to be within 1-2%. 
This is is significantly lower than any other data that can 



be found in the literature (e.g. Kanada & Aoyama (19831). 
As can be seen, variation of an input parameter (e.g. s or 
0i„) changes both, SNR and v. Eq. Q suggests that SNR is 
composed of visibility and a coherent term, therefore deter- 
mining the direct S2Vf?-dependence is complicated. By sep- 
arate inspection of these two terms we can measure v, but 



a: 

to 



40 - 



50 



?0 



10 



0E 



+ 



f , = 300mm 
f, = 800mm 
model 



0.0 0.1 0.2 0.3 0.4 0.5 0.6 

P 2 

Figure 9. Signal-to-noise times visibility vs. p 2 for fi = 300 mm 
and 800 mm data (error bars are in the same range for both data 
sets, but for clarity reasons only represented for the fx = 800 mm- 
case). The theoretical prediction corresponds to M = 900 modes 
as can be obtained from an independent speckle-count measure- 
ment (see fig. 1 1 [ ) . This is consistent with the six data points at 
p 2 :C0.2. Deviations occur for SNR > 1 200 (dotted line) as is the 
case for p 2 = 0.25 and 0.5, discussion see text. 



SNRcoh is difficult to measure and not directly ascertainable 
from the experimental data. However, rather than measur- 
ing each of its components separately, we divide the task of 
determining a concise SNR-mode\ into two steps: 

(i) Verifying the coherent model: We start out with the 
assumption that SNR is factorisable. In order to show that 
this assumption is true, it is sufficient to show that v x SNR 
gives the same response as if illuminating the spectrograph 
with monochromatic light (from eq. [3] SNR co h = v x SNR, 
this is the factorisation). 

(ii) Determine visibility dependence: After showing 



that (i) is fulfilled, the problem of modal noise prediction 
reduces to the problem of predicting the visibility. The vis- 
ibility needs to be investigated. 



4.1 Verifying the coherent model 

We show that the experimental results support that SNR 

v x SNR can be identified 



is factorisable as stated in (i) 
with SNR co h because the product follows the the functional 
dependence as described by eq. |l]| for (a) p 2 and (b) the 
mode number M. 

Fig.jsjjshows (v x SNR) vs p 2 (~ aperture area). The co- 
herent model is supported by the six data points at p 2 ^ 0.2 
of the high resolution (/i = 800 mm) and the low resolution 
(/i = 300 mm) setup. The theoretical prediction assumes a 
mode volume M of 900 modes. Not included into the theo- 
retical fit is the p 2 = 0.25 and 0.5 data points as their low 
value indicates a systematic deviation from the theory. We 
therefore expect an additional source of noise that becomes 
effective when approaching high SNR -values (of the order 
of 1 200 and higher, see dotted line in fig. |9| and longer 
exposure times (see also fig. [5] where the SNR decreases af- 
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Figure 10. Number of guided modes vs. &i n . M^f is derived from 
modal noise and visibility measurements as in eq. [9] The dotted 
line represents the theoretical fit to <x 2 of the two central values 
at 0i n =0.08 and 0.13. Counting speckles under laser illumination, 
the ^independence obeys the square law (solid line) but with an 
offset. 



ter reaching a maximum) - For an investigation with lower 
.S/Vi?- values, higher p 2 -values, a higher resolution setup or a 
test-fibre guiding fewer modes would be required. The error 
bars in fig. [9] are dominated by the uncertainty in visibility 
prediction. 

The p 2 -values are derived from the different slit sizes re- 
stricting the beam at the spectrograph entrance (compare 
fig-0: 

2 = Ad ^ sD = 4s_ = 4s 
P ~ A i ~ (tt/4)D 2 ~ tyD ~ nfa/Fout 1 ' 



where in the second step the approximation Ad ~ sD has 
been used. With slit width s, diameter of the far field image 
D, and fo =4 mm (focal length of the lens imaging the fibre 
far-field onto the spectrograph entrance slit). The input fo- 
cal ratio is Fi n = 3.9. Due to FRD the corresponding output 
focal ratio F out is 3.5. With eq. Q the available slit-widths 
0.1, 0.2, 0.25 and 0.5 mm thus correspond to p 2 = 0.1, 0.2, 
0.25 and 0.5, respectively (p 2 =0.05 and 0.15 are realized by 
restricting the 0.1 mm and the 0.2 mm slit in cross-dispersion 
direction, see fig. |4(b)[ ). The other input/output focal ratios 
are treated similarly. 

In further investigation for the additional source of noise, 
we measured the SNR for the no-slit case, this could not 
fully explain the additional noise that causes the p 2 = 0.5 
data points to deviate from the predicted value. This in turn 
might indicate that this instability arises from the speckle 
pattern itself possibly being subject to changes for long in- 
tegration times. A slight decline could be detected for longer 
integration periods and using a neutral density filter (20 % 
transmission) accordingly. Especially in the very sensitive 
high STVR-regime the impact can become very significant. 

To further approve the assumption in (i) i.e. 
SNRaoh^v x SNR we investigate its ^-dependence. If SNR 
is factorisable as proposed in eq. |3| and if we believe in the 
p 2 -dependence as stated; using eq. |l]) the value SNR X v is 



a measure for the number of modes that are excited in the 
fibre at the central wavelength: 



M eff = (SNR ■ 



1-P' 



(9) 



This is a function of p 2 (known value, set by experi- 
ment) and the SNR x u-product. Hence it follows Meff = 
M eff (SNR x v). 

Fig. [10] shows effective number of excited modes M e // vs. 
angle of the input cone 9i n , where the crosses are the M e //- 
values (see eq. |9|) which are derived from SNR x v-values. 
The dashed line is a fit on the experimental data oc 9 2 n and 
the diamonds represent the number of modes derived from 
laser illumination. The fit on the speckle data (solid line) 
obeys the square law, but has an offset. This implies that 
there will always be a certain number of modes, irrespective 
how small we choose the input angle 9i n - 
Note here that p 2 varies with the size of the far field at the 
output of the fibre, which is subject to input focal ratio and 
focal ratio degradation. Moreover, the term effective mode 
number is used here to comply with the fact that the modes 
are not equally excited. The total number of excited modes 
is lower than the number of modes that is allowed for the 
mode volume in eq. |2|. Depending on the symmetry of the 
fibre launch, different mode groups are excited. 



According to Daino et al. (19801, the number of excited 
modes is equal to the number of speckles that can be 
counted when observing the fibre under monochromatic 
light. A standard He-Ne laser instead of the white light 
source. Here it has to be considered that the laser light 
is polarized, thus the expected number of speckles is only 
half compared to the case of unpolarized input light. The 
number of modes depicted in fig. [10] is therefore derived by 
multiplying the number of counted speckles by two. 
The relationship between input focal ratios f/3.9 and 
f/6.5 (6in = 0.13 and 0.08) is conform with the theory, 
i.e. the number of excited modes rises with the square of 
the coupling angle. However, the 8i n = 0.04 and 0.19 are 
inconsistent with the theory. However, all M e //-values show 
consistence with the number of modes as derived from 
counting the speckles under laser illumination, indicating 
that the product SNR x v can indeed be identified with 
the number of guided modes and thus supporting the 



factorization described in (i) 

It is not straightforward to decide whether to base the 
model prediction on input focal ratio, output focal ratio 
or a mixture of both, because part of the focal ratio 
degradation could occur at the each of the fibre end- 



faces (compare length-effect, Poppett & Allington-Smith 
( |2010 1 ) , leading to a complex mode-population depen- 
dence along the fibre. We chose modal noise to be a 
function of 9% n as this seems to fit the coherent model 
best. For a better grounded statement, the actual mode 
population in the fibre as a function of the coupling angle 
needs to be known; Unfortunately, this is hard to determine. 
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Figure 11. Visibility values for different input focal ratios. The 
experimental values well reflect a (9 _2 -law (dashed line). 
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Figure 12. Visibility values for different slit widths, compared 
for two different resolving powers (/i = 300 mm and 800 mm). 



4.2 Determining visibility dependence 

The goal of this section is to determine the visibility depen- 
dence of the spectrograph parameters; input focal ratio and 
slit- width at the spectrograph entrance. Previous consider- 
ations have shown that the initial hypothesis holds true; 
the visibility signal-to-noise product is consistent with the 
coherent model as described in eq. (TTj) . Fu rther more, the 
number of excited modes rises with 8 (fig. 10 1 as deter- 



mined from directly counting the number of speckles at the 
fibre exit (due to the offset slightly different from the the- 
oretical prediction). Therefore we can conclude that, once 
the visibility is known, the modal noise is determined (see 
section 4.1 1. Thus, the visibility is now subject to closer in- 



vestigation. A fibre that guides a large number of modes will 
consequently possess a large range of propagation constants 
and thus will exhibit a low visibility. Likewise, a wide slit 
reduces visibility (cf. eq. |6]l). 

The graph in fig. [II] shows the visibility as a function of 
the input focal ratio 9in (same data as used in fig. [8] and 
fig. 10 1. From eq. |6| we expect the visibility to follow a 



0~ -law, which indeed reflects the experimental data if we 
accept an offset of 1.7%. This offset becomes dominant 
for fast input beams (8 > 0.05 rad). Visibility as function 
of slit width is depicted in fig. |12| Although eq. |6| pre- 
dicts an s~ 1 -dependence, the experimental values suggest 
s — ' 5 -dependence . These first results give us a parameter- 
ization, which reflects notably well the predictions of the 
simplified model. However, it also shows that the analytical 
model needs further refinement and may best be understood 
by simulation based investigations that reflect the complex 
physical processes when speckle images of adjacent wave- 
lengths interfere at the detector plane. 

Counter-intuitively, the supposed higher spectrograph 
resolving power for fx = 800 mm (= higher coherence) 
shows lower visibility values than the lower resolving power 
setup (/i = 300 mm). It cannot be excluded that this effect 
is an inherent feature of the way the spectra overlap at 
the detector. However, the authors hold the view that the 
apparent independence is more likely due to problems in 



the setup: repeated measurements with re-alignment of the 
fibre-launch resulted in a SWR-degradation of about 15%; 
also differences in the optical performance of the applied 
collimator lenses possibly affect the measurement. In 
order to conclude for the visibility dependence on spectral 
resolving power this point needs clarification. 



5 CONCLUSIONS AND OUTLOOK 

Modal noise is a phenomenon of increasing importance as 
it poses a fundamental constraint on the achievable signal- 
to-noise of future high resolution spectroscopy applications. 
In this paper, special focus has been drawn on the char- 
acteristics of fibre-spectrographs used in astronomy, which 
we found to be insufficiently addressed in the literature. Fur- 
thermore, we experimentally verify that the model for modal 
noise prediction can be separated out into a visibility and a 
coherent term. Thus, modal noise prediction reduces mainly 
to the prediction of the visibility. 

We have outlined the background understanding of modal 
noise as presented in existing literature, along with a sim- 
ple ray-propagation model for visibility prediction. A test- 
bench for experimental investigations was built, with resolv- 
ing powers in the range of R — 10 000 — 200 000. Its design 
and data-analysis are presented. The results support the 
separation approach, showing that modal noise factorizes 
as predicted. As the coherent term can be determined from 
speckle-statistic considerations, it is sufficient to know the 
visibility for noise prediction. Therefore, our further exami- 
nations concentrated on predicting the visibility. The func- 
tional dependence of the experimentally derived visibility 
values show partial agreement with the simplified theoretical 
description. These considerations need yet to be extended by 
further simulation work, while carefully paying attention to 
the spectrograph geometries. This rather extensive simula- 
tion work will potentially complete the modal noise picture 
which will be examined in a future paper. 
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fs = V X 



V2 



An estimation for /„ is derived as follows; 
v c _ c 
R RXc 

Where v c and A c are the central frequency and the central 
wavelength of the spectral resolution element with width 
Av = v\ — vi. Hence, 

L6 2 



fs/fc ~ 2n c RX 
can be substituted in eq. Q, yielding 



1 + 2 



2n c RX 



-1/2 



Note: It has to be said that skew rays are omitted here for 
simplification. 



APPENDIX B: MODAL NOISE FOR 
PARTIALLY POLARIZED LIGHT 

For a circular fibre the number of excited modes is equal in 
both polarization planes: 

M 



Mm = Mi = 



(Bl) 



The two speckle patterns for each polarization plane inter- 
fere independently and are integrated additively at the de- 
tector. We can ascribe two independent signal-to-noise val- 
ues to each of these patterns. 



£BVR||,x = a - a/mJ 



where a = pj(v^f\- 



(B2) 



If the intensities of these two 



planes differ (partial linear polarization) by 
I\\-I±=p-I (B3) 
(where p is polarization degree); the combined SNR results 



APPENDIX A: VISIBILITY ESTIMATION 

The bandwidth f c in eq. [4] can be estimated as follows: As- 
suming a light cone of focal ratio F = 1/26 is coupled into 
the fibre, the difference in propagation constant between on- 
axis rays (Po) and rays coupled with the angle 6 (propaga- 
tion constant /3) is given by 



Af3 = f3 -p oy /l-eyn 2 c 



2 raj 



•#> = 



2nu 



2 Re c/n c 



where n c is the core refractive index, c the speed of light in 
vacuum and v is the frequency of the electromagnetic wave, 
Po = 2w/(c/n c ) has been used. 

The time needed for a mode traveling the whole fibre length 
L is given by t = (3L/2kv. Accordingly, the time delay 
between on-axis and off-axis rays (launched at 9) is then 
At = A/3L/2-KV. Identifying this time delay with f^T yields 



fc 



2-kv 



2n c c 



A/3L L6 2 



SNR 



SNR 2 + (l-p)SNR 2 L = 
y /(M i + (l-p)M ± ) = ayfl 



(B4) 



In the last step eq. |B1| has been used. For 10 % polar- 
iz ed light, the measured SNR will be reduced by 1 - 
yjl - 0.1/2 ^ 2 % compared to non-polarized light. 
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